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ABSTRACT 

We investigate colour selection techniques for high redshift galaxies in the UKIDSS Ul- 
tra Deep Survey Early Data Release (UDS EDR). Combined with very deep Subaru optical 
photometry, the depth (K AB = 22.5) and area (0.62deg 2 ) of the UDS EDR allows us to in- 
vestigate optical/near-IR selection using a large sample of over 30,000 objects. By using the 
B — z, z — K colour-colour diagram (the BzK technique) we identify over 7500 candidate 
galaxies at z > 1.4, which can be further separated into passive and starforming systems 
(pBzK and sBzK respectively). Our unique sample allows us to identify a new feature not 
previously seen in BzK diagrams, consistent with the passively evolving track of early type 
galaxies at z < 1.4. We also compare the BzK technique with the R — K colour selection of 
Extremely Red Objects (EROs) and the J — K selection of Distant Red Galaxies (DRGs), and 
quantify the overlap between these populations. We find that the majority of DRGs, at these 
relatively bright magnitudes are also EROs. Since previous studies have found that DRGs at 
these magnitudes have redshifts of z ~ 1 we determine that these DRG/ERO galaxies have 
SEDs consistent with being dusty star- forming galaxies or AGN at z < 2. Finally we observe 
a flattening in the number counts of pBzK galaxies, similar to other studies, which may in- 
dicate that we are sampling the luminosity function of passive z > 1 galaxies over a narrow 
redshift range. 
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1 INTRODUCTION 

Colour selection can provide an efficient technique to identify high- 
redshift galaxies based on known rest-frame SED features. How- 
ever, most of these colour criteria are sensitive to more than one 
type of SED feature, i.e. more than one galaxy type and at poten- 
tially different redshifts. 

EROs are one class of colour- selected galaxies, these are con- 
sistent with either dust-reddened star forming galaxies or passive 
galaxies at z > 1 (e.g. |Cimatti et aT]|2002| |Roche et af||2002| 
Simpson et al. 2006]). Extending these techniques to higher red- 
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shifts Da ddi et al.| f2004) suggested a more refined criteria based 
on two colour cuts in B — z and z — K, now known as the BzK 
technique, where BzK = (z! — K)ab — (B — z')ab- Star form- 
ing galaxies in the range 1.4 < z < 2.5 are selected by requiring 
BzK ^ —0.2. This was determined by analysing the BzK prop- 
erties of z > 1.4 star-forming galaxies, classified through [Oil] 
emission and UV spectra. Passive galaxies in the same redshift 
range require BzK < —0.2 and [z — K)ab > 2.5. This was 
determined by analysing the BzK properties of spectroscopically 
confirmed passive galaxies at z > 1.4. Star-forming and passive 
BzK galaxies will now be referred to as sBzK and pBzK respec- 
tively. Finally the selection of galaxies using near-infrared colours 
(J — K)ab > 1.3 (DRGs) has been shown to be an effective tech- 
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nique for selecting galaxies out to z > 2 based on the B aimer break 
( |Franx et al.|2003] >. 

The effectiveness of these various photometric techniques, 
however, is still unqualified, due in part to the small sample sizes, 
or shallow depths, used in previous follow-up studies. The latest 
generation of large area, deep near-infrared surveys, such as the 
UKIDSS UDS, now provide a representative sample of colour- 
selected sources. These can be used to determine the proportions 
of galaxy type selected using different colour criteria, their redshift 
ranges and the overlaps between different techniques. 

Throughout this paper we use a concordance cosmology with 
Q m 0.3, Q A 0.7, and H = 70 km s" 1 Mpc -1 . The mag- 
nitudes and colours quoted above are based on 2 arcsec diame- 
ter aperture magnitudes, as are all magnitudes and colours quoted 
throughout this paper. 



2 COLOUR SELECTION OF HIGH-Z GALAXIES 
2.1 Data Set and Sample Definitions 

The UKIRT Infrared Deep Sky Survey (UKIDSS) has been running 
since spring 2005 and comprises 5 sub- surveys covering different 
areas and depths ( Lawrence et al. 2006 ). The UKIDSS survey uses 
the Wide-Field Camera (Casali et al., in preparation) on the United 
Kingdom Infrared Telescope (UKIRT). This study makes extensive 
use of J and if-band data from the Ultra Deep Survey early data 
release (the UDS EDR; see |Dye et al.|2006| ). The UDS is the deep- 
est of the 5 UKIDSS sub-surveys and aims to reach a final depth 
of K A b = 25.0, H AB = 25.4, J A b = 26.0 (5a, point source) 
over an area of 0.8 deg 2 . The size of the UDS field significantly 
reduces the effects of cosmic variance and on this scale the UDS 
will be the deepest image ever produced, providing an unparalleled 
number of candidate high-redshift sources. Due to a small change 
in the UDS field centre shortly after the beginning of the survey, the 
current JK imaging is not uniform over the entire 0.8 deg 2 field. 
Consequently, for the purposes of this paper, we restrict our analy- 
sis to the 0.6 deg 2 central region which has uniform JK data and 
reaches depths of Kab — 22.5 and Jab — 22.5 (5cr limits). By 
conducting simulations, the completeness at this magnitude limit 
is determined to be above 70% for point sources. For details of 
the completeness estimation, image stacking, mosaicing and cat- 
alog extraction procedures see |Foucaud et al.| ( |2006| >. In addition 
to J and if -band data from the UDS, deep B (B Urn = 28.2), 
V (Vum = 27.5), R (Rum = 27.6), i (i Um = 27.5) and z '- 
band (z' Urn — 26.5) data from Subaru Suprimecam are also used 
( |Sekiguchi et al.|2005| ). 

The wealth of multi-band data available in this field, together 
with their depth and large area coverage, has enabled the construc- 
tion of a highly-detailed BzK colour-colour plot which displays 
many interesting features including a new branch (Figure p}. Fo r 
consistency the standard BzK definitions of Dad di et al.| (2004), 
as discussed in the introduction, were used to construct our BzK 
sample. K band sources (> 5cr) were crossmatched with our op- 
tical Subaru i band source catalog; since the Subaru field is not 
entirely coincident with the UDS field this reduced the usable area. 
For K band sources (> 5cr) that did not have an optical match, 2 
arcsec aperture magnitudes were obtained directly from the Subaru 
images. All BzK sources were then visually verified to make sure 
they were not caused by diffraction spikes from saturated sources or 
cross-talk effects ( |Dye et al. 2006 ). Additionally, only sources that 
were not within the halos of saturated optical sources were used, 




Figure 1. Plot showing the BzK colour-colour diagram for sources in the 
UDS EDR. The sBzK and pBzK regions are marked accordingly. Branch 1 
is the new branch of galaxies found in this study and Branch 2 is the stellar 
branch. Overlaid are no evolution model tracks for a range of SED types, 
starting at z=0 and ending at z=1.5. 



which left an area of 0.5591 deg . All if-matched sources were 
used for BzK selection unless both the B and z band magnitudes 
were < la limits (0.7% of the sample, outside of saturated ha- 
los), since these cannot be constrained within the BzK plane. This 
resulted in the selection of 6736 sBzK and 816 pBzK. 

In addition to BzK selection the availability of R and J band 
data in this field has also enabled the extraction of EROs and DRGs. 
In this case the ERO selection was defined as (R — K) vega > 5.3 
and is carried out using the prescription presented in Simpson et al. 
pOOoT ), resulting in 4621 EROs (to K A b = 22.5). ERO selection 
was carried out independently because Simps onet al.| ([2006) uses 
a different method of optical-NIR matching than used here. The 
DRG sample used here is constructed in Foucaud et al. (2006 ) using 
a selection criteria of (J — K)ab > 1.3, resulting in 369 DRGs 
(to Kab = 21.2). Only those DRGs within the Subaru region and 
outside of saturated halos are used, this leaves 330 DRGs suitable 
for cross matching with NIR-optical sources. 



2.2 BzK Number Counts 

The differential if -band number counts for the different popula- 
tions in our survey (Figure [2] Table [TJ show that star-forming BzKs 
increase in number sharply toward fainter magnitudes. In contrast 
pBzKs seem to exhibit a knee at K ~ 21 where the number counts 
clearly turn over before the limiting if -band magnitude, this fea- 
ture is also seen in |Kong et al.| ( [200 6). The small redshift range 
over which pBzKs are thought to be selected (Kong et al. 2006) 
could provide an explanation for the knee and the relatively faint 
magnitude at which pBzK number counts begin. That is, that the 
number counts are directly sampling the luminosity function of this 
population. 

We note that the galaxies forming this knee feature have very 
faint B-band magnitudes, close to the completeness limit in this 
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Figure 2. K-band differential number counts for all sources in the field 
and BzK sources. The sBzK population rises rapidly to our completeness 
limit, while the pBzK counts exhibit a knee at Kab ~ 21. Data from the 
Deep3a-F sample of Kong et al. ( 2006) is included for comparison. 

band (Bu m = 28.2). However fainter B magnitudes will push the 
B — z colours to the red and hence these sources should not be 
missed by our pBzK selection. At Kab ~ 21 completeness is 
^100% ( [Foucaud et al.|2 006) so the feature is unlikely to be due 
to incompleteness. We conclude that this knee feature is likely to be 
real. We note that the turn-over corresponds to absolute magnitudes 
Mr ~ —23.6 at z = 1.4, which is close to the value of M* 
determined for the if -band luminosity function at these redshifts 
( |Cirasuolo et al.|2006) . 

As can be seen in Figure[2]our results are similar in form to the 
data from the Deep3a-F survey ( Kong et al. 2006), especially for 
the overall source number counts. However, for the sBzK number 
counts there is a slight magnitude offset, most likely due to cosmic 
variance. The UDS EDR has greater dynamic range in K due to our 
greater depth and area. This dynamic range will increase toward 
even fainter magnitudes with further UDS releases. 

2.3 BzK New Branch Galaxies 

Due to the large number of detected sources in this survey field 
(34098 in total), and the comparative depth of the optical data avail- 
able, there are a large number of sources available for BzK analysis. 
It is clear from Figure [TJthat a new feature is visible running paral- 
lel with the stellar branch, at the bottom of the plot, but with redder 
z — K colour, roughly defined by the region: 

0.3(5 - z') - 0.2 < — K) < 0.3(B - z') + 0.4 (1) 

with clear branch separation at (B — z') > 2.5. 

Also plotted on Figure [T] are no-evolution model tracks for 
different types of galaxy. These were created by redshifting SEDs 
from King & Ellis (1985), convolved with the appropriate filter re- 
sponse curve to get the relevant band magnitudes from which the 
(B—z) and (z'—K) colours can be calculated. The predicted track 
for early type galaxies (E/S0) corresponds very well with the new 
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Table 1. Differential Number Counts in log(N/deg 2 /0.5mag) bins for sBzK 
and pBzK in the UDS EDR. 
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branch of galaxies, especially at low redshifts. As can be seen, at 
higher redshift galaxies are no longer populating the predicted E/S0 
model track. This is partly due to incompleteness at these high red- 
shifts but also because these model tracks do not take into account 
evolution, so early type galaxies will be bluer at higher redshifts 
than the models predict. As such they only provide a guide as to 
where different types of galaxies exist within the BzK plane. The 
tracks of later-type spirals fall predominantly within the main trunk 
of galaxies in the BzK diagram and all lie within the sBzK region 
by z^l.5. The proximity of the Sab model track to that of E/S0 
suggests that there may be some contamination of the early type 
branch by spirals, especially at low redshift. 



3 ERO/BZK COMPARISON 

The large number of EROs, BzKs and DRGs provided by the UDS 
EDR for the first time allows the overlaps between these popula- 
tions to be analysed in detail. Previous studies have been made of 
the overlaps between EROs and BzKs, using smaller numbers of 
sources than are available here. Kong et al. ( 2006) find that 41% of 
EROs are also selected as BzKs in their Deep3a-F data, and 29% 
from their Daddi-F data. This compares with Da ddi et al.| (2004) 
who find that ^35% of their ERO selected sources are also BzK se- 
lected. In our much larger study we find that 60.6 ± 1.5% of EROs 
(to Kab = 22.5) are selected as BzK, of either variety, which is a 
higher fraction than found in previous studies. This is not surprising 
since the UDS EDR is deeper and contains fainter EROs, of which 
a larger fraction are likely to represent z > 1.4 objects than would 
be expected in a brighter sample. 

Figure [3] shows the overlaps between the four populations 
studied here and the number of sources involved in each case. To 
ensure we are only comparing samples to the same completeness 
depth only EROs and BzKs selected to Kab = 21.2 are used in 
this diagram, since this is the limit for DRG selection. However, 
when comparing EROs and BzKs we also quote numbers to a depth 
of Kab = 22.5 in parentheses in the following discussion. As can 
be seen, 204(792) out of 214(816) pBzKs are additionally selected 
as EROs. The reason for this strong overlap is clear from the posi- 
tion of EROs on the BzK diagram (Figure[4]); EROs are located in a 
clearly defined region with red colours in z' — K and B — z' which 
covers the entire region occupied by pBzKs. It is also clear from 
Figure [4] that sBzKs are far less preferentially selected as EROs - 
283(2007) out of 892(6736) sBzKs, or 31.7 ± 1.4%(29.8 ± 1.2%). 



© RAS, MNRAS 000,[T}{5] 



4 K. Lane et al. 



sBzK (892) pBzK (214) 





DRG (330) 
49 


577 1 


10 




32 ■ 
-1— 

85 1 







98 


66 


j 712 




198 1 


138 



ERO (1297) 



Figure 3. Venn diagram showing the overlaps in the number of objects clas- 
sified as EROs, BzKs and DRGs. All numbers given here are to the depth 
of the DRG sample (K AB = 21.2). 

From the position of BzKs on a plot of R — K colour versus K 
magnitude (Figure[5]) it can be seen that there is a clear difference in 
distribution between pBzKs and sBzKs. pBzKs have redder R — K 
colour than sBzKs, with a narrow range around R — K — 4.6 
(FWHM~0.24, see side panel in Figure [5}. This is in contrast to 
sBzKs which tend to be broad in R — K colour but occupy regions 
with fainter K magnitude (see top panel in Figure[5]). This could be 
explained if sBzKs constitute the higher redshift (z > 1.4) end of 
the ERO population. Using the photo-z catalog presented in Cira- 
|suolo et a l. (2006 ) shows the same redshift behaviour. However, 
given that these photo-z values are based on the same photometric 
data used here to construct the BzK population, and use the same 
SED features used to place galaxies at z ~ 1.4 this does not essen- 
tially add any new information, but serves as a consistency check. 



4 DRG/ERO COMPARISON 

The DRG selection technique was originally intended to preferen- 
tially select galaxies at z > 2. However the most recent generation 
of large area NIR surveys (e.g. Grazian et al. 2006 Conselice et al. 
2006) are starting to probe the bright, low redshift, end of the DRG 
population. Based on photometric redshifts Quadri et al. ( 2007} find 
> 70% of DRGs (to K vega < 20) from the Multiwavelength Sur- 
vey by Yale-Chile (MUSYC) to lie at z > 1.8. Compared to the 
UDS EDR, however, the greater depth and substantially smaller 
area used in Quadri et al. (2007) means that a much higher fraction 
of the DRGs selected are likely to lie at higher redshifts. A more 
comparable study is Conselice et al. ( 2006) in which a DRG sam- 
ple is produced to a depth of K vega = 20.5 over a slightly larger 
area than used in this study. They find their DRG sample to peak 
at z rsj 1.2, when using photometric redshifts, and at z ~ 1.0, 
when using spectroscopic redshifts, with only 4% at z > 2. This 
redshift distribution does not change when their DRG sample is 
cut to the same depth as the DRG sample constructed in this study 
(Kvega = 19.3). We therefore take z ~ 1.0 to also be the likely 




(B-z')ab 

Figure 4. Figure showing the BzK colours of EROs. Almost all pBzKs are 
EROs as well as a large fraction of sBzKs. This can also be seen in Figure 

® 

redshift distribution of our DRG sample. A new picture is starting 
to emerge in which the DRG criterion not only selects z > 2 pas- 
sive galaxies, but also galaxies at redshifts of z ~ 1. It is therefore 
interesting to see how DRGs cross match with photometric popu- 
lations of a more established nature, such as EROs and BzKs. 

Most DRGs in the UDS EDR are also selected as EROs 
(75.5 =b 6.3%, see Figure [3}. This means that these galaxies are 
red in both R — K and J — K colours which implies they have a 
steep slope in their SEDs across this colour range. An SED feature 
of this kind, at z ~ 1, is indicative of dusty star-forming galaxies 
or AGN. XMM-Newton X-ray data is available for this field (Ueda 
et al. in prep.), but only 6 (1.8 ± 0.7%) of our DRGs are found to 
be X-ray detected. At the depth of the XMM data, however, this 
is only sufficient to rule out luminous AGN with relatively modest 
Nh columns (e.g. Nh < 10 23 cm -2 , assuming Lx — 10 44 erg 
s _1 at z = 1.5). Therefore the presence of highly obscured and/or 
lower luminosity AGN cannot be ruled out. Of the joint ERO and 
DRG galaxies, 34.1 ± 4.3% are sBzK and 26.5 ± 3.7% are pBzK. 
The number of pBzKs in this overlap group suggests that a size- 
able portion of these pBzKs are likely to be either obscured AGN 
or dusty star forming galaxies rather than purely "passive" pBzK. 

It would therefore appear that DRGs not only select z > 2 
galaxies but also form part of the z~l ERO population, based 
in part on the DRG redshift determination from |Conselice et al.| 
( 2006), as dust reddened star forming galaxies as well as a fraction 
of possible obscured AGN (see also Cons elice et al. 2006). 



5 DRG/BZK COMPARISON 

Approximately half of the galaxies selected as DRGs are also se- 
lected as BzKs (117 sBzK and 66 pBzK). Although DRGs are se- 
lected as BzKs across the full DRG K magnitude and J — K colour 
range, they tend to have fainter K magnitudes and are more effi- 
ciently selected by redder J — K colours. Of these joint selections, 
72.6 ± 10.4% of the sBzK selected DRGs were also selected as 
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Figure 5. An R — K diagram showing BzK selected sources that are also 
selected as EROs. The histogram to the right shows the R — K colour of 
EROs and ERO selected p and sBzK. It is clear that pBzKs have redder 
colours, whereas sBzKs cover a broad range in R — K colour. The top 
histogram shows X-band aperture magnitudes for the same subsets. Much 
like the overall ERO population, sBzKs tend to be at higher magnitudes and 
peak around Kab ~ 22, conversely pBzKs have a much broader distribu- 
tion which peaks at Kab ~ 21.2. 



EROs, as were all of the pBzK selected DRGs. This is as expected, 
due to the large fraction of DRGs selected as EROs (section[4). 

|Reddy et aL|p005] > find that -30% of DRGs are found to be 
sBzKs, which is similar to the fraction of DRGs found as sBzKs 
in our study. The study also finds that —10% of sBzKs are DRGs, 
as are —30% of pBzKs; though sample errors will be high since 
only 17 pBzKs are used. This is in good agreement with this paper 
which finds that 13.1 ± 1.3% of sBzKs and 30.8 ± 4.3% of pBzKs 
are also selected as DRGs. The large area of the UDS EDR data 
is probing the lower redshift, brighter end of the DRG luminosity 
function, so a large overlap with BzK- selected galaxies would be 
expected (we observe —55%). 

The larger fraction of pBzKs selected as DRGs could be due 
to the narrower and lower redshift range of pBzKs, as hypothesised 
in previous studies, combined with the DRG technique selecting 
galaxies at z — 1 . It could also be due to some shared astrophysical 
feature, such as AGN (as discussed in section]?}. Whichever is cor- 
rect, it is likely that BzKs constitute at least the z — 1.4 — 2.0 part 
of the DRG population. Those DRG not selected as BzK are likely 
to be at z < 1.4, as in Daddi et al.| ( [2004| ), though some could be at 
z > 2.5. 



6 CONCLUSIONS 

For the first time a statistically significant study of the overlaps 
between ERO, BzK and DRG galaxy populations has been carried 
out. We compare 1297 EROs, 330 DRGs and 1106 BzKs, selected 
from 0.559 1 deg 2 of imaging to K A b = 21.2. It is found that BzKs 
are consistent with being the z > 1.4 end of the ERO population, as 
would be expected from the definition of these selection techniques. 
It is becoming clear from the new generation of large area 



surveys that the DRG selection criterion is not only effective at 
extracting high z galaxies, but also dusty star forming galaxies and 
obscured AGN, at z—1, particularly at brighter if-band magnitudes 
(Kab < 22). In this study we find that 249 of 330 DRGs are also 
selected as EROs. Those DRGs also selected as pBzKs tend to be at 
the red end of the DRG sample while those also selected as sBzKs 
have a range in J — K colour but are faint in if -band magnitude. 
This is consistent with DRGs selecting galaxies over a broad red- 
shift range, from EROs at z — 1 to BzKs at z — 1.5. Deeper UDS 
data will also allow us to probe the faint end of the DRG selection 
regime to ascertain the effectiveness of this technique at selecting 
galaxies at z > 2. 

The depth of our Subaru optical data has allowed us to deter- 
mine the number-magnitude relations of the BzK samples. sBzKs 
have a broad range in magnitude, consistent with a broad range 
in redshift and/or luminosity. The pBzKs, however, exhibit a turn 
over in their number counts that are consistent with pBzKs inhab- 
iting a narrow redshift range at z — 1.5. The UDS EDR combined 
with this Subaru data has also allowed the construction of the most 
complete BzK diagram yet seen. The number of sources has al- 
lowed the identification of a new feature that is most likely to be 
the no-evolution track of passive early type galaxies with increas- 
ing redshift. 
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